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Quantum signatures of classical multifractal measures 

Moritz Sclionwetter 1 and Eduardo G. Altmann 1 
1 Max Planck Institute for the Physics of Complex Systems, 01187 Dresden, Germany 

A clear signature of classical chaoticity in the quantum regime is the fractal Weyl law, which 
connects the density of eigenstates to the dimension Do of the classical invariant set of open systems. 
Quantum systems of interest are often partially open (e.g., cavities in which trajectories are partially 
reflected/absorbed). In the corresponding classical systems Do is trivial (equal to the phase-space 
dimension), and the fractality is manifested in the (multifractal) spectrum of Renyi dimensions 
D q . In this paper we investigate the effect of such multifractality on the Weyl law. Our numerical 
simulations in area-preserving maps show for a wide range of configurations and system sizes M that 
(i) the Weyl law is governed by a dimension different from D 0 = 2 and (ii) the observed dimension 
oscillates as a function of M and other relevant parameters. We propose a classical model which 
considers an undersampled measure of the chaotic invariant set, explains our two observations, 
and predicts that the Weyl law is governed by a non-trivial dimension D asymp totic < Do in the 
semi-classical limit M —> oo. 

PACS numbers: 05.45.Mt,05.45.-a,05.45.Df 


I. INTRODUCTION 


Fractality is a well-known signature of chaotic dynam¬ 
ics m- In open quantum systems the fractality of the 
classical dynamics appears in the distribution of eigen¬ 
states [4] and in the Weyl law [5]. Recent studies consid¬ 
ered the effect of periodic orbits [1] and weak chaos CHS 
and were conducted on 4-dimensional Hamiltonian sys¬ 
tems j7] and different area-preserving maps [0 [8HTT)1 . In 
these systems the fractal Weyl law can be written as 

N(\Vi\ > Hcutoff) ~ M D o S> /2, (1) 


es) 

where M is the dimension of the Hilbert space, D’ is 
the fractal dimension of the chaotic saddle, and |iz| C utoff 
is a cutoff separating long-lived from short-lived states 
with eigenvalues v t . 

Optical cavities m and other relevant physical sys¬ 
tems Il2j are not completely open but instead show par¬ 
tial absorption. A simple physical picture is a billiard 
in which trajectories hitting the boundary are partially 
reflected and partially absorbed. In contrast to open sys¬ 
tems (full absorption), in partially absorbing systems, the 
box-counting dimension D 0 is an integer number equal to 
the phase-space dimension m Fractality in the classi¬ 
cal dynamics of such systems is still present mi and can 
be quantified using the spectrum of Renyi dimensions [1] 
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where p ( x ) is the relevant probability measure 
(f n dp ( x ) = 1) in the phase-space Q and the sum is over 
e-sized boxes Bi (with 0 < pt := f B p(x)dx) used to 
partition the phase-space. For q = 0, D 0 is recovered. A 
distribution p ( x ) with a non-trivial dependence of D q on 
q is called multifractal. 

In this paper we investigate the quantum mechanical 
signatures of multifractality in the classical phase-space. 


While a straightforward application of the fractal Weyl 
law to partially open systems predicts a trivial scaling 
(N ~ M), as also argued in [Tl] . non-trivial scalings have 
been reported in numerical investigations m In order 
to investigate the robustness of fractality in the Weyl law 
we study it at the transition between monofractal and 
multifractal. This is done by varying the reflectivity R > 
0 of an absorbing region localized in the phase-space. For 
all nonzero reflectivities we obtain non-trivial Weyl laws 
with effective dimensions that show strong oscillations as 
a function of system size and other parameters. 

We argue that these observations are due to an effective 
undersampling of the classical measure in the quantum 
regime and that the non-trivial scaling remains valid in 
the semiclassical limit M —> oo. 


II. FRACTALS IN CLASSICAL DYNAMICS 

As numerically convenient examples of Hamiltonian 
chaotic systems we consider here fully chaotic symplectic 
maps of the form /:i „4 x n +i, defined in a bounded 
phase-space x G $7. To incorporate the effect of partial 
absorption we associate to each trajectory an intensity 
In, with I 0 = 1, I n+ 1 = R(x n )I n , and R(x) G [0,1] is a 
real valued function describing the reflectivity of different 
regions in the phase-space. 

In order to investigate the differences between partially 
and fully open systems it is instructive to consider the 
special case of R(x) = R < 1 in a small absorbing region 
of the phase-space (i.e. R(x) = 1 everywhere else). If 
R = 0 in the absorbing region the system corresponds to 
usual open maps and the set of initial conditions xq with 
liuin-^oo I n = 1 form the forward trapped set of the map 
(the stable manifold of the chaotic saddle [2] |T5] [IB]). In 
this case the measure is uniformly distributed on a fractal 
support and the Renyi spectrum is flat ( D q = Do). In 
analogy, for a partially absorbing system 0 < R < 1 the 
measure in a phase-space region Bi (e.g., the i- th box of 
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a) baker map b) cat map 
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FIG. 1: (Color online) Classical dynamics of area-preserving maps ( x n ,y n ) (x n +i, Vn+i) with absorbing regions. The pairs 
of plots a) and b) show the action of one application of the baker’s map (x n +i, 2 /„+i) = (3a; — [3a;,,J, y n /3 + \3y n \/3) and the 
cat map (xn+i, Vn+i) = (2x„ + y n mod l,3x n + 2y n mod 1), respectively. The areas surrounded by dotted lines in the left 
plots indicate the absorbing regions (R < 1) of each map. The dark regions in the right plots show the first images of the 
absorbing region. 



a partition by a regular grid) should be proportional to 
the average intensity /„ for n —> oo of randomly drawn 
xq £ Bi. The finite time normalised measure of Bi after 
n iterations can thus be written as 

_ (I)n,i /o', 

/ T\ ’ 

Z-^iX 1 /n,i 

where the average intensity (...)n,i is computed over ini¬ 
tial conditions xq £ Bi. This definition is consistent with 
the usual definition of the stable manifold measure pQ 
considering the intensities I £ [0,1] of trajectories as 
weights (see also Ref. ED- 

As representative cases of fully chaotic area-preserving 
maps we investigate in detail the baker Q] and the cat 
map nu as defined in Fig. [I] In the case of the baker 
map the absorbing region was chosen in order to allow 
for comparison with analytical calculations. In the case 
of the cat map we chose an arbitrary asymmetric region. 
The essential point for our investigation of the role of 
partial opening is that both choices lead to a non-trivial 
Dq for R = 0. In Fig. [2] we report the distribution of the 
measure p in the phase-space for a particular reflectivity 
and the spectrum D q . For the baker map the calcu¬ 
lations can be done analytically (see Appendix 0 and 
confirm the validity of our numerical methods of esti¬ 
mating [ii and D q . In both maps Dq = 2 and D q is 
non-trivial, confirming the result of Ref. [13] that the 
apparent self-similar distribution of p observed in the 
phase-space is properly quantified only through the full 
multifractal spectrum D q . 


III. FRACTALS IN QUANTUM SYSTEMS 

In this section we introduce the quantised version of 
the maps of Sec. |TT] and compute Weyl’s law (0- By 
open quantum map we mean a propagator m 

M = UP, (4) 



FIG. 2: (Color online) Multifractal properties of the maps 
defined in Fig. [T] Top row: finite time approximation of a 
typical phase-space distribution of the measure /r, defined in 
Eq. j3j. Bottom row: Renyi spectrum D q for R = 0.3 (see 
Appendix[A]for the analytical calculations on the baker map). 


where IA is the quantised closed system and V is a pro¬ 
jector that contains all the information about absorption. 
In the case of full absorption, V removes the part of the 
quantum wave-function ip{x) that overlaps with the ab¬ 
sorbing region, but does not alter the remaining part. In 
the case of partial absorption, V reduces the intensity bv 


a factor R(x) as 


(x\V\ 


= R[x) \i/j(x)\ 2 . See App. 


B 


for the matrix representations of the maps used in this 
paper. 


In an M-dimensional Hilbert space (effective H = ) 

maps A4 have have M left (and right) eigenstates with 
complex energies e.; = £) + iF»/2. A4 acts like a time 
evolution operator so its eigenvalues Vi are related to 
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a) quantum baker map b) quantum cat map 
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FIG. 3: Distribution of eigenvalues \v\, see Eq. |5|)in the 
quantised baker (a) and cat (b) maps. Two matrix sizes were 
used: M = 3' (gray line) and M = 3 9 (black line) [T§]. The 
vertical dashed line marks the cutoff value used for the fractal 
Weyl law calculations in Fig. [djid. 
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Vi = exp 

r~n e \ 

= exp 
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Since r, £ [Rq and £) £ IR it follows that M| £ [0,1]. 
Hereafter we use \v l \ to quantify the lifetime of the state, 
e.g. long-lived states have M| 1$ 1 (I\ <C 1). 

We calculated the eigenvalues of A4 for both the baker 
and cat map for different reflectivities R and matrix sizes 
M = 3 n with n £ IN [Till . The normalised distributions of 
M| are shown in Fig. [3] The states concentrate around 
a typical value M typ , which is M typ ~ 0 for R = 0> 
increases with R , and can be approximated as the median 
of p(M). For increasing matrix-sizes M the distribution 
becomes more concentrated around v | typ . To investigate 
this tendency in further detail, and in line with the fractal 
Weyl law, we concentrate on the number of eigenstates 
with \Ui\ > M cutoff . 

Figures [4jrd show the relative number of eigenstates 
N(yi : \vi\ > |i/| cuto fi)/M as a function of the matrix- 
size M. For this plot, Mcutoff was chosen in such a way 
that Mcutoff > Mt yp but still small enough to allow the 
calculation of the scaling for a wide range of M, see ver¬ 
tical dashed lines in Fig. 1E3- The division by M allows 


for an improved analysis of the scaling, which according 
to the Weyl law prediction ([I]) should be do — 1, where 
do = D 0 — 1 is the partial fractal dimension of the forward 
trapped set along the unstable direction. D 0 denotes the 
full fractal dimension of this set. It is connected to the 
dimension of the chaotic saddle D^ by [21 

D (S) 

D^ S) =2D 0 ~2 O *, = -§-■ (6) 

The usual fractal Weyl law is confirmed in Fig. [4jtd for 
the case R = 0. For small R , the curves follow the R = 0 
case for small M before deviating from it. With increas¬ 
ing R, the deviation is observed already for smaller M. 
Interestingly, even for large R and M the numerical re¬ 
sults do not approach a constant value as predicted by 
the classical dimension do = 1. Instead, a non-trivial 
scaling with an effective dimension between 0 and 1 is 
observed. Even if the numerical results do not allow for a 
conclusion regarding the scaling in the limit M —> oo, the 
roughly constant scaling for a broad range of values in M 
shows the importance of the effective dimension (similar 
effective dimensions have been employed to investigate 
classical [20j and quantum [8] Hamiltonian systems with 
mixed phase-space). 

In order to more carefully analyse the scaling we de¬ 
fine the local dimension di oc as the slope of the N/M 
vs. M curves between two nearby matrix sizes M. Fig¬ 
ures [TJae show the dependency of d\ oc on the choice of 
the cutoff M cutoff . For small cutoffs Mcutoff < M typ 
we observe di oc = 1, in agreement with the observa¬ 
tion in Fig. [3] that the distribution of eigenvalues be¬ 
comes concentrated around M typ (almost all states count 
as long-lived). A more interesting behavior appears for 
Mcutoff > Mt yp > w bere at least two regimes governed 
by non-trivial di oc can be identified before di oc starts to 
fluctuate for Mcutoff ► 1 due to a lack of statistics: (II) 
an regime of intermediate M cutoff with strong (almost 
periodic) oscillations of d\ oc ; and (12) a regime in which 
dioc ~ dQ R ~°\ These regimes are particularly visible for 
small R (because v \ typ —I 0), but can be identified also 
for larger R. Surprisingly, Figs. S f show that regime 
II increases for increasing M (strong oscillations cover a 
wider range of cutoffs). This suggests that the oscillating 
(local) scaling is not an artifact of small M but instead 
that it will prevail in the large M limit. 

In summary, our numerical observations of the quan¬ 
tum maps show surprising features. In particular: (i) 
the classical prediction do = 1 was not observed; (ii) the 
effective scaling varies smoothly with R; and (iii) strong 
oscillations we observed in di oc vs. M cu toff- 


IV. SAMPLING MULTI-FRACTAL MEASURES 

In this section we propose an explanation for the scal¬ 
ing of the number of long-lived states with M observed 
in the previous section based on the classical dynamics 
of the system. 
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FIG. 4: (Color online) Scaling behaviour of the quantum states for the baker (a-c) and cat (d-f) maps. Top row: Fraction 
of eigenvalues above the cutoff value M cutofI as a function of the system size M at different reflectivities. The dotted lines 
mark the expected scaling for R = 0 (dg R ^ 0 ^ = ln(2)/ln(3) for the baker map and d < ' 0 R ^°' > ss 0.42 for the cat map). Middle 
row: Local scaling exponent di oc (M) := (lniV(M) — lnA r (M/3)) /ln3, where N(M) denotes N(vi : \vi\ > M cutoff ), at M = 3 9 
for different values of I? as a function of M cutoff (colors and symbols as in top row). We show results up to a cutoff where 
N(M ) < 0.005M after which point the lack of resonances leads to strong oscillations of di oc . The vertical dashed lines indicate 
the cutoff values used in the top row. The small ticks on the upper axis mark the positions of the medians of the eigenvalue 
distributions (ss |z/| t ). Bottom row: d\ oc {M) as a function of the cutoff for M = {3 7 ,3 8 ,3 9 } (lighter shades correspond to 
smaller M). 
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A. General Argument 

We first consider how the usual fractal Weyl law is 
obtained from the classical phase-space of fully absorb¬ 
ing maps M- Lu et al. [5] argue that the states can 
be considered as non-overlapping and that each of them 
covers an area of h = 1/M. This localisation happens 
in the forward trapped set (for the left eigenstates |55J) 
mm, which is smooth. Therefore, the number of dif¬ 
ferent long-lived states depends on the distribution of 
the forward trapped set in the perpendicular direction 
(along the backward trapped set or, equivalently, along 
the unstable manifold of the chaotic saddle). The rea¬ 
soning above leads to the following constructive classical 
procedure: 

1. Cover an intersection of the forward trapped set 
with segments of size h = e = 1/M. 

2. Count the number A/ of boxes (segments) needed 
to cover the forward trapped set, which estimates 
the phase-space volume available to the long-lived 
eigenstates. 

3. The number of long-lived states N(M ) is propor¬ 
tional to Nf,. 

Since A/, ~ M d ° we obtain N(M) ~ M d °, the fractal 
Weyl’s law 0. The scaling exponent do is the par- 
tial fractal dimension of the forward trapped set along 
the unstable direction <©■ Schomerus and Tworzydlo 0 
provide an alternative explanation based on the phase- 
space areas that do not escape up to a given time (which 
grows with M). Both explanations are equivalent be¬ 
cause asymptotically the A-sized boxes intersecting the 
invariant set correspond to the non-escaping phase-space 
areas. 

We now generalize the ideas above to the case of par¬ 
tial absorption. In this case we expect for large M the 
long-lived states to be distributed according to the clas¬ 
sical measure p of the forward trapped set. Therefore, 
it is essential to consider a procedure that accounts for 
the non-uniformity of p. In line with the usual interpre¬ 
tation of p as a probability measure, we modify point 2. 
above and attribute to each segment i a probability Pi 
of being counted as long-lived (i.e., a probability that a 
long-lived state is localized in it). This probability has to 
be large (small) where the measure /q of the segment i is 
large (small) because states localise on high measure ar¬ 
eas of the phase-space. At the same time, for large M the 
number of possibly long-lived states grows and pi should 
grow because of the multiple chances of one state being 
localized in segment i. The simplest stochastic process 
agreeing with these two intuitions corresponds to per¬ 
forming S independent trials with a success-probability 
Pi- This can be incorporated in the procedure described 
above by replacing point 2. by 

2’. Count each segment i as part of the ac¬ 
cessible phase-space volume with a probability 


Pi = 1 — (1 — Pi) s , where (1 — pi) s is the proba¬ 
bility that the box is rejected S times. 

For one realization of this process the available volume 
V ( M) is proportional to the number Nj, of counted seg¬ 
ments. Therefore, point 3. is generalised to 3’. as 

M M 

AI(M)<xE[iV 6 ] = ^p i = ^l-(l- / x i ) s , (7) 

i —1 i=1 

where E[...] is the expected value computed over mul¬ 
tiple sampling realisations and the number of trials S 
is associated to the total number of states S = SqM. 
In the case of a fully absorbing region, p t = 1/M 
for all M boxes intersecting the forward trapped set 
and pi = 0 everywhere else. Eq. ( 0 ) simplifies to 

E [Nb] = M (l — (1 — 1/M) S ) S ~ y °°> J\f anc l we recover 
the fractal Weyl law. Analogously, in the fully reflecting 

case Pi = 1/M and E[JVj] S ~ ¥oa > M and we recover the 
traditional Weyl law. 

At the heart of the reasoning above is a crucial dif¬ 
ference between the computation of do and of the Weyl 
scaling, do is obtained by first taking the limit of sam¬ 
pled trajectories S to infinity S —> oo (construction of 
the measure) and afterwards computing the scaling of 
filled boxes in the limit of small box sizes e —► 0. In 
the Weyl law the analogues of these two limits (S —> oo 
and £ —>• 0) happen simultaneously when M —> oo be¬ 
cause the maximum number of eigenstates grows with A1 
( S oc M) and the phase-space resolution increases with 
M (h = e = 1/M). Therefore, the scaling in the Weyl 
law corresponds to the estimation of do when sampling 
the measure with a sample size 

S = S 0 M = S'oE -1 (8) 

in ( 0 . The proportionality constant Sq effectively con¬ 
trols how many segments (long-lived states) will be found 
for the different M values and therefore plays a similar 
role as the cutoff value M cu toff i n the quantum case. 


B. Application to Baker Map 


We now apply the general sampling method proposed 
in the previous section to the baker map (defined in 
Fig. 0 . The expected number of long-lived boxes E [AT&] 
for £ = 3~ n can be computed analytically (see Appendix 

[0 


E [N b \ 


3 n -J2 2, 


-i 


1=0 


l\(n — l)\ 



R l \ 
(2 + R) n J 


s 

. (9) 


In Fig. 0 we plot Eq. (0) and its local derivatives as 
a function of its three parameters £, So, and R. The 
similarities to the quantum observations shown in Fig. 0 
are remarkable (the comparison is based on the identi¬ 
fications £ = h = 1/M , N{vi : \vi\ > M cuto ff) oc E [N b ], 
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FIG. 5: (Color online) Scaling behaviour of the classical 
expected values for the baker map. The black dotted line 
denotes cIq 101 . a) Rescaled expected value of accepted boxes 
at an arbitrary sample size of S(M) = SoM with So = 0.3. 
E[1V(,] grows at fixed e monotonically with R. The shaded 
areas cover the range between E[A},] ± a[Nb], where u is the 
standard deviation, b) Scaling exponent d\ oc (M = 3 9 ) as a 
function of 1 /So. The lower the value of R, the lower the 
curve intersects the dashed vertical line (the So choice used 
in a), c) dioc(M = 3 n ) with n € {7,8,9}. Lighter shades 
correspond to smaller M. 



FIG. 6: (Color online) Comparison of the quantum data to 
the classical model. The symbols show the number of long- 
lived states. The solid lines are the classical values of aE[IV{,] 
obtained by choosing So and a for each R so that the right¬ 
most two points of the quantum data are matched as close as 
possible. 


and the similar role played by So and Mcutoff' In partic¬ 
ular, we reproduce the three main observations reported 
in Sec. |III| (i) the scaling of the curves is different from 
do = 1; (ii) for increasing R, the effective scaling devi¬ 
ates more and more from and approaches a non¬ 

trivial dimension (Fig. [5^,); (iii) the local slope di oc show 
strong oscillations with So (Fig. [5 }d) in a range of So 
values which increases with M (Fig. l5t) . One quantum 
feature which is not immediately apparent in Fig. is 
the plateau of di oc ~ d , the regime 12 mentioned at 
the end of Sec. m To understand this we need to take a 
closer look at the mapping from So to the cutoff Mcutoff - 
The limiting case So —► oo can be safely identified with 
| v | cutoff —> 0. In this limit our model correctly reproduces 
the quantum observations (di oc H► 1). The identification 
of the opposing limit So —> 0 with M cutoff —>• 1 is also 
possible but less meaningful, since in this limit both the 
number of resonances, and 7V b become very small, and as 
a result the determination of di oc becomes prone to fluc¬ 
tuations. Between these limits, we can assume that 1/So 
grows monotonically with Mcutoff' However, the func¬ 
tional form relating 1/So and M cu toff unknown (e.g., 
it depends on the distribution of eigenstates) which does 
not allow for a careful quantitative comparison of the two 
plots (e.g., a non-linear transformation of the x-axis in 

Fig. [HJo could generate the plateaus around d^ - ^ seen 
in Figs. [djod. 

We now search for a quantitative comparison between 
our classical sampling model and the quantum observa¬ 
tions. The only parameter we have to fit is So in the 
sampling model, which according to the arguments of 
Sec . |IV A| should reproduce the scaling of the long-lived 
states with large M. Accordingly, for each value of R 
we choose the value of So which reproduces the local 









































7 


slope between the last two N/M points in the quantum 
data (Figs. |4^,d) and we shift the E [iV&] curve obtained 
to match the rightmost point of each quantum data¬ 
set (this last shift fixes the proportionality factor a in 
N(M ) = aE[iV(,]). We restrict So to the range [3 _1 ,3], 
so that S is of the order of M since these values give the 
best agreement with the whole range of quantum data. 
Because of the oscillations in d\ oc there may be more 
than one choice of So [ 21 ] • The results for the different 
R shown in Fig. [6] confirm that our model describes the 
data also for values of M much smaller than the region 
used in fixing So, a. 

The quantitative success of our sampling model mo¬ 
tivates us to consider its behavior in the semiclassical 
limit M —> oo. In Fig. 0 we plot Eq. ([9]) in an increased 
range of M = 1/e values using for the different R curves 
the value of So used in Fig.JhJo (estimated from the quan¬ 
tum data). We again see the presence of log-periodic os¬ 
cillations (see Ref. [22] for a review on log-periodicities 
in fractals). While these oscillations play a major role 
in dioc and for the range of M accessible in the quan¬ 
tum computations, they do not hinder the estimation of 
an asymptotic dimension dasymptotic as the general trend 
of the curves. For R = 0 the results confirm that for 
dasymptotic = d^^ and that our sampling correctly ac¬ 
counts for the usual fractal Weyl law. More interestingly, 
dasymptotic < 1 for all R, strongly suggesting that the 
non-trivial scaling di oc observed in the quantum data are 
not finite M effects but persist in the asymptotic limit. 
This implies that the fractal Weyl law of partially open 
systems does not follow the classical prediction do = 1. 
Comparing the R dependence of d aS y mp totic with the an¬ 
alytical calculated d q indicates that roughly q ss 0.3 can 
be used to estimate this scaling. 


V. DISCUSSIONS 

In summary, we have investigated the statistical prop¬ 
erties of long-lived eigenstates of chaotic systems with 
partial absorption. Contrary to the naive prediction ob¬ 
tained combining the classical fractal dimension do and 
the fractal Weyl law, we observe numerically that for 
a broad range of parameters the number of long-lived 
states grows sub-linearly with the dimensionality M of 
the Hilbert space. We explain these observations by 
considering a model which is consistent with the quan¬ 
tum to classical correspondence in the semiclassical limit 
M —> oo but that accounts for an undersampling of the 
classical measure at finite M. This undersampling com¬ 
bined with the strong spatial fluctuations of the multi¬ 
fractal classical measure is responsible for the effective 
reduction of the available phase-space volume for long- 
lived eigenstates. 

Our most surprising finding is that even the asymp¬ 
totic (M —► oo) scaling of the number of long-lived states 
differ from the do = 1 prediction. This result is compat¬ 
ible with our numerical simulations at least for a broad 



M 



FIG. 7: (Color online) Asymptotic scaling of the 

classical prediction. a) Extrapolation of E [Nb] of the 
curves in Fig. [6^ without the shift (i.e. a = 1). R 
increases from the black (bottom) to the brown (top) 
curve. The dashed lines follow the asymptotic scaling 

M d asymptotic- 1 wit h daaymptotic = C™ ^2-hiM^ ^ where 

the average is computed over all possible combinations 
{(Ad, M') : M' < M A (Ad, M') G {3 10 ,3 11 ,..., 3 40 } 2 }. In b) 
we show these dimensions together with their standard de¬ 
viations. In grey are analytically obtained d q curves, see 
Eq. (IA13 1. Note, that for fixed R and q G IR + , do > d q > doo 

ID- 


range of cutoffs in \ v\. In any practical situation the limit 
M —> oc is not achievable and the relevance of our model 
is that it allows to understand the effective dimension ob¬ 
served quantum mechanically. We have also considered 
alternative models in which we assume that the eigen¬ 
states concentrate exclusively on the phase-space regions 
with highest measure /x, or on boxes with a /x-dependent 
probability that we compute without including a sam¬ 
pling process. None of these models is able to reproduce 
the quantum observations, see Appendix [D] 

Our results are not limited to area-preserving maps 
or to systems with localized absorbing regions. Instead, 
they should be visible in systems of any dimension for 
which a multifractal spectrum of a classical measure /i 
appears. In n-dimensional systems the long-lived reso¬ 
nances are expected to concentrate in n-cuboids of vol¬ 
ume h and the sampling has to be carried out on a hy- 























perplane intersecting these volumes. In general partially- 
absorbing systems there is no analogous fully open sys¬ 
tem (R 0 for the cases above). In these cases, the 
scaling of Weyl’s law is expected to follow the dimen¬ 
sion of the regions of highest probability [3D] for small M 
(analogous to do of the R = 0 case in the examples above) 
and show a transition towards the dimension d asymp t 0 tic 
of the undersampled fi (for the large M case). The multi- 
fractality in the localization of wave functions can appear 
also without a clear classical counterpart, e.g. in the case 
of phase transitions in disordered systems )23| [Mj ■ It 
would be interesting to investigate whether the number 
of such states show similar dependencies on system size 
as the ones we observed. 
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Appendix A: Analytical D q in the baker map 

Assume a ternary baker map as defined in Fig JT] with 
reflectivities Rl, in the left, Rm in the middle and in 
the right vertical strip. In this system it is possible to 
calculate the number of necessary iterations to get an ex¬ 
act result up to some minimum £ m ; n . The size of features 
(stripes) in the ternary baker map after n iterations is at 
least 1/3". If £ m ; n is smaller than that value then the 
box-counting will pick up the features as areas and D 0 
will approach 2. 

Let £ [0,1]. After the first iteration the inten¬ 

sity in the top horizontal third is Rl,, the central strip 
has an intensity of Rm, and the bottom third carries Ar . 
The total intensity is 

Itot.i = Rl + Rm + Rr- (Al) 

The next iteration gives the following total intensity: 

hot,2 = Rl, + RlRm + RlRr + 

RmRl + Rm + RmRr + (A2) 

RrRl + RrRm + R'n 

and so forth. At the n-th step we obtain the total inten¬ 
sity 

n n—k / \ 

w = (A3) 

= (Rl + Rm + Rr) 

Here we used the fact that the number of strips with 
a certain combination of reflectivities R^R l M R^ k ^ 1 is 
given by the multinomial coefficient (//) := k\i\(n-k-i)\ • 

Now we focus on the multifractal dimension spectrum 
([2]) of the distribution generated by this map. First, since 


the measure is constant in the horizontal direction, we 
can restrict ourselves to a vertical line and calculate the 
partial dimensions d q along this line. Assume this line 
is covered by disjoint segments of length e = 3 _ ". The 
dimensions of the full distributions are then simply ob¬ 
tained from 


with 


D q = 1 


lim 


ln Ei=i M? 


q 1 — q n-t 00 n ln 3 


(A4) 

(A5) 


Here ^ is the intensity on the i- th segment normalised 
by hot,n, so that 


& i = 1 


2=1 

Or, again using multinomial factors: 

n n—k / 

Tl 


with 


k =0 1=0 


t^n;k,l — 


U i)^ k ' 1 


= 1 


p/e pZ TDn—k—l 
-Miot ,n 


With that we can write the sum in (A5) as 


(A6) 


(A7) 


(AS) 


1 -!>?= 

(A9) 

n n—k ( n \ ( p/e pZ p n—k—l\Q 

1 \k,l) ) 

n hh (Rl+rm+rrT* 

(A10) 

(Rl + Rl + R^ n 

(Rl + Rm + RR) nq 

(All) 

1 ( r l + + Rr) 

(Rl + Rm + RrY 

(A12) 


Thus we obtain for d a 


d q = - lim 

1 1 — q n—too 


1 (^l+^m+^r) 

n 111 (Rl+Rm+Rr ) 9 

n ln 3 


(A13) 


This is independent of n so we can drop the limit and get 
ln m±K)_ 


dq — 


(.Rl+-Rm+-Rr ) 9 


(1 — g) In 3 ' (AW) 

The limit q —> 1 can be obtained using L’Hopital’s rule. 

1 InEflxM? 


di := lim 
=4- d -\ = 


g->i 1 — q n ln 3 
In (Rl + Rm + Rr) 
ln3 


(A15) 


Rl In Rl + Rm In Rm + Rr In Rr 
(Rl + Rm + Rr) ln/3) 
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Asymptotically d q goes to 


,. , In (Rl + RM + Rr) ~ In (Umax) 

Inn d q = - — -. 

9->oo 1 In d 

where R max := max{/? L , R M , Rr}- 


(A16) 


Appendix B: Quantised maps 


The matrix representation of the propagator of the 
quantised version of the baker map defined in Fig. [l] is 


baker 


/Ibaker 


fl-M/3 0 0 

0 \/R1m/3 0 

o 0 1 M / 3. 


(Bl) 


with 


/Ibaker b 


M 


Fm/3 0 0 

0 R M /3 o 

0 0 F m / 3 


(B2) 


where F M . m ,n = ^ exp [-^(n + l/2)(?n + 1/2)] are 
Fourier transformations from position to momentum rep¬ 
resentation and F// denotes the inverse transformation. 
For an explanation of a quantisation scheme see e.g. S3] . 

Analogously, the quantised cat map has the following 
form: 


A^cat — 


M 


ljkf/6 0 0 

o VR^-m /2 o 

o o l M /3_ 


(B3) 


where M is even and the entries are 




t/M 


exp 


/ o 2\ 

—— (m — mn + n 
M v ' 


(B4) 


Appendix C: Derivation of E[N b 


Starting from the general equation |7| for the expected 
number of accepted boxes here we show how to arrive at 
the formula for the baker map ®>- 

The i-th box is found with a probability Pi = 1 — (1 — 
Pi) s . Analogous to (A8) we can restrict ourselves to 
different boxes and define p n -.k,i := 1 — (1 — p n -k,i) S ■ The 
expected value of the number of found boxes is 


n n—k 


=Eft=xm 


k =0 1=0 


kj 


Pn;k,l • 


(Cl) 


Plugging in the expression for p n - y k,i we get 

E [N b ] = (C2) 

V V ( n \(^ - (^ - R £ Rl M R K k - 1 2 
V (Rl + Rm + Rr)" 



Equation follows from setting Rl,r = 1, Rm = R, 
carrying out the inner sum, and renaming the indices. 

Appendix D: Discussion of alternatives to sampling 


Here we want to show that taking into account only 
the regions of phase-space with the highest measures by 
a cutoff in the measure can not reproduce the scaling 
observed in the quantum maps. 

First, let us consider a fixed cutoff. In this case it is 
easy to see that for any value of p cu tofi = const, there 
exists a n* with 


Tn;k,l ^ Mcutoff: ^ tl • 

Therefore do n ^°° > o. 

The second case is an n-dependent cutoff value 
Mcutoff,n = j° u )° ff with constant I cu toff- In this case if 
all Rl,m,r < 1 there also exists a n* with p n -k,i < 
ftcutoff.nj Vn > n*. In other words: do "~ > ' 00 > 0. An ex¬ 
ception here arises if Rl,r = 1 and Rm £ [0,1). Then, 
for all 0 < /cutoff < 1 after a number of iterations n* 
all middle thirds fall below the cutoff do In 2/In 3 

from above. 

So, by associating the number of long-lived eigenstates 
of boxes above a given threshold we obtain an effective di¬ 
mension contradicting our numerical observations. This 
indicates that the (large) regions with low measure p still 
contain substantial portions of the eigenstates. 

The third alternative we want to discuss is the case 
where each box contributes to the long-living phase-space 
volume proportionally to its intensity 


Pn-k : l CC In;k,l — Pn;k,lRot,n- (El) 

As before we identify M = 3 n . In this context we need to 
limit the contribution of each box. p n -k t i = minjl, al n -k,i} 
and compute E [N b \ = Jfk=o J2?=o (k,i) Pn:k >i‘ The factor 
a plays a similar role as the cutoff in the sense that for 
a —> oo (a —i 0) we count all (zero) boxes. This method 
reproduces the correct scalings in the cases Rl,m,r = 1, 
and Rl,r = 1,Rm = 0. However, it fails to reproduce the 
increase of the effective dimension with R and M and the 
oscillating behaviour of di oc with Hcutoff- 
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